A Consistent Picture for Resonance Neutron Peak and ARPES Spectra in High-T c 
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The spectra observed in ARPES measurements are exam- 
ined together with the resonance peak observed in neutron 
scatterings, based on the slave-boson approach to t — t' — J 
model. We show that the peak/dip/hump features arise 
from the scattering of electrons by collective spin excitations 
which, at the same time, give rise to the neutron resonance 
mode. The doping dependences and the dispersions of the 
peak/dip/hump positions are shown to be consistent with ex- 
periments. Our results indicate that recently observed cos(66 l ) 
deviation from the pure d-wave also result from the renormal- 
ization by spin fluctuations. 

PACS number: 71.10.-w,74.25.Jb,74.25.Ha,79.60.Fy 



Both angle resolved photoemission spectroscopy 
(ARPES) and neutron scattering experiments have 
played important roles in the studies of high-T c supercon- 
ductors. It has been shown by ARPES that the spectral 
lineshape possesses a peak/dip/hump structure in the su- 
perconducting(SC) state ^J,^|. The anomalous momen- 
tum, temperature and doping dependences of the spec- 
tral lineshape suggest that electrons are strongly cou- 
pled to collective excitations centered at (tt, tt) and these 
collective excitations are related to the pairing interac- 
tion ||. On the other hand, the most prominent feature 
of the spin susceptibility observed in neutron scattering 
studies is the sharp resonance peak at (tt, tt) in the su- 
perconducting or pseudogap states (j ||. It has been 
speculated HQ, from a comparison of these two kinds 
of experimental data, that the collective excitations are 
the resonance modes in the neutron scattering experi- 
ments. This idea has been further explored qualitatively 
in a phenomenological spin-fermion model ||, in which 
the resonance mode is identified as the propagating col- 
lective spin excitations and the scattering of electrons 
by these spin modes gives rise to the anomalous spectral 
lineshape. Using the slave-boson theory for the t — t 1 — J 
model, Brinckmann and Lee Jy] have investigated the spin 
resonance and its evolution with doping. 

Recent ARPES study reveals the specular dependences 
of the peak, dip and hump energies with doping JlC| ] . Fur- 
thermore, a small cos(6(9) (9 is the Fermi surface angle) 
deviation from the pure d-wave structure is observed [ |ll| . 
These results and their correlation with the mode energy 



inferred from neutron data are the essential ingredient 
for a consistent picture of the ARPES lineshape and the 
neutron data. In this paper, we examine these issues and 
show that they can be quantitatively reproduced based 
on the slave-boson approach to the 2D t — t' — J model. 
The model reads, 



<IJ>,<7 



<ij>',tr 



(1) 



where < ij > denotes the nearest- neighbor (n.n.) bond 
and < ij >' the next n.n. bond. In the slave-boson 
method, the physical electron operators Ci a are expressed 
by slave bosons bi carrying the charge and fermions fi a 
representing the spin; Ci a = bffi a . We consider the 
d-wave SC state with the order parameters Ay =< 
fafol ~ fnfjl > and Xij = E CT < ftJi° >■ in which 
bosons condense 6, —>< bi >= \/5 (6 is the hole concen- 
tration) |12j. Then, the mean-field Hamiltonian of Eq.(l) 



H m = J2 ZkftJk* - Y, M/fcV-fcl + cc) + 2NJ'(xl + Ag), 
k 



ka 
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where the dispersion for fermions = —2(St + 
J'xo)[ cos (kx) + cos(k y )] — 4St' cos(k x ) cos(k y ) — /i, and the 
gap A(fc) = 2J'A [cos(fc x ) - cos(k y )], with J' = 3J/8. 
The mean-field parameters xoi and the chemical po- 
tential fj, for different doping 5 are obtained from a self- 
consistent calculation p2| . 

In the above mean-field calculation, we have set < Si > 
=0 p|. In order to consider the response to external 
magnetic and electric fields, we include the antiferromag- 
netic (AF) fluctuation by writing H = H m + H', where 
H' = H — H m , and treating H' as a perturbation to 
H m . In the Hartree-Fock approximation, this reproduces 
the mean-field results fPj|| . To investigate the AF fluc- 
tuations, we consider the scattering of electrons off spin 
fluctuations in H' . As a first step, we calculate the spin 
susceptibility in the random-phase approximation (RPA) 
as shown in Fig. 1(a), 



x(q,w) = xo(q,w)/[l + aJ(q)xo(q,u>)]. 



(3) 
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Here J(q) = J(cosq x + cosq y ) and Xo(q, w ) is the un- 
perturbed spin susceptibility which is calculated from 
the fermionic bubbles representing particle-hole excita- 
tions. Following Ref. ||, we choose a = 0.34 in or- 
der to set the AF instability at S = 0.02, which is 
the experimental observed value. This is the only ad- 
justed parameter throughout the paper. The fermionic 
self-energy is obtained from the lowest-order contribu- 
tion of the scatterings of fermions off spin fluctuations. 
In the SC state, there are two different self-energies 
E s and Su, as shown in Fig. 1(b) and (c), which renor- 
malize the fermionic dispersion and the SC gap, re- 
spectively. In the previous study only E s is in- 
cluded in their calculations. We will show that the 
inclusion of £„, lead the otherwise flat doping depen- 
dence of the ARPES peak to be consistent with experi- 
ments. The fermionic Green's function is calculated by 
G/(k,w) = [G7 x (k,w) + (A fe + S tu ) 2 G7 1 (-k,-^)]-' 1 
with G/ (k,w) = [iuj - e k - E s (k, u))]~ l . In the SC 
state, bosons condense and the physical electron Green's 
function can be approximated by G(k, u>) « <5G/(q, uj) 
. Then, the spectral function of electrons is calculated 
from A(k, to) — — (l/7r)ImG(k, to). Numerical calcula- 
tions are performed at low temperature T = 0.005 J, with 
t = 2 J, t' = -0.45*, and J = 0.13eV. 

We first analyze the imaginary part of the spin sus- 
ceptibilities at Q = (tt, 7r). It develops sharp reso- 
nance peaks for various doping densities. This result 
has been reported in Ref. ||. In the framework of 
the (i-wave BCS theory, the origin of these peaks has 
been discussed [p| Jl^ , [l5| . Essentially the peak arises from 
a collective spin excitation mode corresponding to 1 + 
aJ(Q)Rexo(Q,^) — and negligibly small Imxo(Q,^)- 
It is due to a step-like rise of Imxo at its gap edge and 
then a logarithmic singularity in Rexo via the Kramers- 
Kroenig relation. This singularity shifts downward the 
collective mode energy and leads it to situate in the spin 
gap, so no damping is expected for the mode. The step- 
like rise in Im^o arises from the flat band(extended van 
Hove singularity) which is observed near (tt, 0) jjll and 
the property that Ak+Q = — Ak for transition momen- 
tum Q due to the <i-wave gap symmetry |H| . 

The lineshape A(q, co)f(ui) (f(u>) is the Fermi distri- 
bution function) of electrons coming from the scatter- 
ings by spin fluctuations for doping S = 0.12 are shown 
in Fig. 2. The results are plotted for several wavevec- 
tors k = (7r, 0), (tt, O.IStt), down to (tt,0.5tt). Clear 
peak/dip/hump structures are present at and near (w, 0). 
In order to understand the origin of the peak/dip/hump 
structure, we plot the self-energy E s [E w is qualitatively 
similar to £ s ] at the Fermi wavevector kp for doping 
5 = 0.12 in the inset of Fig. 2. The solid line denotes 
its real part E s , while the dashed line its imaginary part 
£ s . The corresponding lineshape at kp is expressed by 
the line with open squares in the main panel. Due to the 



coupling to the spin resonance mode, the whole struc- 
ture of the self-energy E s is very similar to that of the 
unperturbed spin susceptibility xo (for a comparison, see 
Fig. 2 in Ref. p5| ). When frequency is below about 
0.5J, E s is equal to zero. Above 0.5J, a step-like rise 
can be seen and is followed by a decrease. Consequently, 
£p has a peak at about the center of the step-like rise of 
E s . The quasiparticle energy is given by the pole of the 
Green's function, which is the solution lu of the equation 
P(k F ,u) = Re[(^-E s (fc^)) 2 -(A fcF +E M ,(fc F ,w)) 2 ] = 
0. We also show P(kp,oj) (dotted line) in the inset of 
Fig. 2. The lowest binding-energy solution of the pole 
equation P(fc F ,w) = is to = —0.42 J. Meanwhile, the 
damping of this mode which is proportional to E s and 
T, w approaches to zero. Therefore, it gives rise to a quasi- 
particle mode which is denoted by the sharp low binding- 
energy peak in the lineshape shown in Fig. 2. As in- 
creases further, the pole equation does not satisfy any- 
more. Near the end of the step-like rise in E s , P(kp,ui) 
reaches its local maximum, meanwhile the imaginary part 
of the self-energy is also near its maximum, therefore 
a dip appears. It indicates that the dip is caused by 
the step-like rise in the imaginary part of the self en- 
ergy. After the dip, P(kp,Lu) decreases with \to\ and near 
to = —1.0 J, it reaches a local minimum. This leads to 
the broad higher binding-energy hump. 

We now discuss the dispersion of the hump/dip/peak 
structure. At the momentum range below and slightly 
above the Fermi wavevector (jr, 0.1571"), the low binding- 
energy peak nearly does not disperse with k. Above 
(tt, 0.37r), the peak starts to move to higher binding- 
energy but cannot move further than that dip seen 
around kp. Then, it disappears gradually with the fur- 
ther increasing of k, due to the unavailable quasiparticle 
states above the Fermi surface at low temperatures. As a 
result, the dip also disappears gradually. These behaviors 
were observed in experiments |fio| . The position of the 
peak determines the normalized gap size due to spin fluc- 
tuations. In Fig. 3, we show a typical dependence of the 
renormalized gap magnitude on the Fermi surface angle 
9, which is the angle between kp and k x . In general, the 
next order deviation from the d-wave starts from cos (69) 
and can be described by the form A(9) = Ai cos(2(9) + 
A2cos(60). However, our results indicate that the best 
fit is A(0) = A max [Scos(26») + (1 - B)cos(60)] with B 
being around 0.88-0.94. This particular form was also 
used in Ref. to fit their ARPES data with roughly 
the same range of B. 

In Fig. 4(a), we show the doping dependence of the po- 
sitions of the hump and the peak, and compare our re- 
sults with the experiment |[o]]. One can see that our 
result for the peak is in reasonably good agreement with 
the experiment. We note that, if we just include the self- 
energy E s in our calculation as was done in Ref. ||, a 
doping independence for the peak positions is obtained. 
Hence the renormalization of the gap by including Y, w is 
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important. The magnitude of the position of the hump 
and its trend with the doping variation are also consistent 
with the experiment. However, the slope of the curve of 
hump versus hole density is flatter than the experimen- 
tal result. Considering that all parameters in our cal- 
culations are chosen according to the well-known values 
and the only adjusted parameter a is fixed by the exper- 
imental observation on AF instability, these results are 
quite satisfactory. The dependence of the ratio between 
the hump and peak energy at (n, 0) on doping concen- 
trations is shown in Fig. 4(b). A flat variation for a wide 
doping level is seen. This result again agrees with the 
experiment 10|. 

An important quantity addressed in the ARPES ex- 
periments is the peak-dip separation, which is shown to 
be close to the mode energy of the neutron peak |jTo| . 
From a comparison of Fig. 2 in Ref. || and Fig. 3, one 
can see that the peak-dip separation in our calculations 
is lower than the neutron peak. However, we would like 
to point out that this separation is sensitive to the ex- 
perimental energy resolution simulated to be T — 0.02J 
here. When T increases, the self-energy E s decreases and 
the quasiparticle peak broads, consequently the peak-dip 
separation increases. However the position of the neu- 
tron scattering peak does not change for different T, it 
just becomes broad. So, this comparison is also sensi- 
tive to r. But, as noted above, the dip stems from the 
step-like edge in E s which is in turn caused by the cou- 
pling to the collective spin mode. On the other hand, this 
spin resonance mode also arises from the step-like edge 
in Irn^o which is due to the extended van Hove singular- 
ity near (jr, 0) and the d-wave symmetry. Therefore, an 
intimate relation between them is suggested. Thus, our 
result seems to provide a consistent picture for the spin 
resonance peak and the hump/dip/peak structure based 
on the spin excitations in a d-wave superconductor. 

Recently, P.J.White et al. |l7| have investigated the 
effect of nonmagnetic Zn impurities on the lineshape 
of Bi2Sr2Ca{Cui^ x Zn x )0$j r $ by ARPES experiment. 
They found that the dip is diminished with Zn doping. 
According to the previous studies by one of the authors 
(Li) jl5) , Zn doping will wash out the van Hove singular- 
ity and cause the decays of quasiparticle states. So, the 
enhancement in Re%o and the step-like rise in Im%o will 
be suppressed. As a result, no clear resonance neutron 
peak appears at certain Zn doping concentration [fl5|| . 
Because the dip is suggested to come from the coupling 
to the spin resonance mode here, the disappearance of 
the dip upon Zn doping may be naturally explained in 
the present framework. 

It is quite encouraging that our results fits various 
kinds of ARPES and neutron scattering data with correct 
trend and reasonable magnitude. There is only one ad- 
justable parameter which is fixed by the experimentally 
observed AF instability. Therefore, our investigation rep- 
resents a natural extension to the work of Brinckmann 



and Lee ||, which addressed the doping dependence of 
the resonance peak in neutron scatterings. In addition, 
in comparison to the qualitative study by Abanov and 
Chubukov Q , our study further furnishes a quantitative 
basis. 

In summary, based on the slave-boson approach to 
the t — t' — J model, we show that the anomalous 
peak/dip/hump structure observed in the ARPES exper- 
iments arises from the coupling of quasiparticle to the col- 
lective spin excitations which gives rise to the resonance 
peak in the neutron scattering experiments. Our inves- 
tigation seems to give a consistent explanation for the 
resonance neutron peak and the ARPES spectra based 
on the spin excitations in a d-wave superconductor. 
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FIGURE CAPTIONS 



Fig.l Feynman diagrams for: (a) the RPA approxima- 
tion to the spin susceptibility coming from particle-hole 
excitations, (b) and (c) the lowest-order contribution to 
the self-energy from fermion-spin excitation scatterings. 

Fig. 2 Photoemission spectra A(k, u>) f (uj) at k = 
(tt, 0), (tt, 0.15tt), (tt, 0.225tt), (tt, 0.3tt), 
(tt, 0.47t), (tt, 0.57r) ( from the up to down lines). Inset 
shows the frequency dependence of the fermion's self- 
energy E s at kp — (7r,0.157r) (Fermi wavevector). The 
solid line corresponds to its real part and the dashed 
line to its imaginary part. Also shown in the inset is 
P(kF, oj)(see text) which is denoted by the dotted line. 

Fig. 3 The renormalized gap size versus the Fermi sur- 
face angle for i5 — 0.08. The open circles are our re- 
sults, while the solid line is the fit using the gap function 
A(0) = A max [Bcos(20) + (1-B) cos(60)] with B = 0.94. 

Fig. 4 Doping dependences of the energy scales for the 
hump and the peak at {%, 0). (a) Doping dependences of 
the hump and peak positions. The open squares and the 
solid triangles are experimental data for the hump and 
peak positions from Ref. |lC(| , respectively, (b) Doping 
dependence of the ratio of the hump position to the peak 
position. 
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Fig. 1 , Li, Mou and Lee 




Fig. 2, Li, Mou and Lee 
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Fig. 3, Li, Mou and Lee 
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